Engineered Optical Nonlocality in Nanostructured Metamaterials 
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We analyze dispersion properties of metal-dielectric nanostructured metamaterials. We demon- 
strate that, in a sharp contrast to the results for the corresponding effective medium, the structure 
demonstrates strong optical nonlocality due to excitation of surface plasmon polaritons that can be 
engineered by changing a ratio between the thicknesses of metal and dielectric layers. In particular, 
this nonlocality allows the existence of an additional extraordinary wave that manifests itself in the 
splitting of the TM-polarized beam scattered at an air-metamaterial interface. 
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It is well accepted that the properties of optical com- 
posites and nanostructured metamaterials can be de- 
scribed by effective parameters derived in the limit when 
structural elements are much smaller than the wave- 
length [1] . The effective medium is an important concept 
of the homogenization theory based on field averaging, 
and it provides a physical insight into the optical response 
of complex micro- and nanostructured media being also 
useful for different types of waves [2]. However, it was 
already established that the effective medium models do 
not provide a complete information, and they should be 
corrected in some cases, e.g. in the recently analyzed case 
of plasmonic nanorod metamaterials in the epsilon-near- 
zero regime, where the performance of such structures is 
affected by nonlocal response [3]. 

Here we demonstrate that the dispersion properties 
of metal-dielectric periodic nanostructured metamateri- 
als are dramatically affected by a nonlocal response due 
to excitation and coupling of surface plasmon polaritons 
at the metal-dielectric interfaces, so in many cases the 
effective medium approach fails to describe correctly the 
optical properties of such structures. The difference is 
dramatic, and it can't be taken into account by small cor- 
rections [3, 4]. Nevertheless, we reveal that the strength 
of optical nonlocality can be engineered in a rather sim- 
ple way, just changing a ratio between the thicknesses 
of metal and dielectric layers, so that the medium with 
equal layers demonstrates weak nonlocality. 

The electromagnetic response of periodic layered 
metal-dielectric nanostructures has been a subject of 
many theoretical and experimental studies. Such struc- 
tures represent the simplest nanostructured metamate- 
rials, they were suggested for a number of applications 
including superlcnscs with subwavelength resolution [5- 
7] , as a simple realization of the so-called hyperlens [8] , as 
well as for nanolithography [9], optical nanocircuits [10], 
invisibility cloaks [11], and even three-dimensional nega- 
tive refraction [12]. In many cases, the effective medium 
is a conventional approach for describing optical prop- 
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FIG. 1: (Color online) Effective medium model. Permittivity 
tensor components as functions of normalized frequency for 
three types of periodic layered metal-dielectric nanostructures 
with the fixed lattice spacing D — di + Gfe, but formed by 
pairs of metal and dielectric layers with the thicknesses (1) 
c?2 = 1.5di, (2) di — c?2, and (3) d\ — 1.5^2, respectively. 



erties of such structures, and it allows to introduce the 
effective permittivities and show that such structures be- 
have as uniaxial metamaterials with the optical axis par- 
allel to the layers. 

In this Letter we reveal that, depending on the 
thicknesses of layers, the metal-dielectric structures can 
demonstrate strong optical nonlocality due to the ex- 
citation of surface plasmon polaritons. In particular, 
we predict the existence of an additional extraordinary 
wave which affect dramatically the scattering of the TM- 
polarized wave at an air-metamaterial interface gener- 
ating two waves with negative and positive refraction, 
respectively. 

We consider a nanostructured metamaterial formed by 
a periodic layered structure (see the insert in Fig. 1), 



2 




kyDIn 



5 6 

Din 



FIG. 2: (Color online) Dispersion diagrams and the geometries of the periodic nanostructures composed of alternating metal 
and dielectric layers. Thicknesses of the layers are: (a) d,2 = 1.5di, (b) di — d2, and (c) d\ = 1.5^2- Two dispersion branches 
corresponding to the actual structures are numbered by Roman numerals. 



where metal and dielectric layers have different thick- 
nesses. When the wavelength of radiation is much larger 
than the thickness of any layer, it is usually assumed that 
the effective medium approach is valid and the permit- 
tivity tensor for an uniaxial anisotropic medium has the 
following form: 
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where £1,62 and d\,d 2 are the dielectric permittivities 
and thicknesses of the layers, respectively. 

A dispersion equation for extraordinary waves sup- 
ported by the such an effective uniaxial medium and 
propagating in the xy-plane has the form 
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and it establishes a relation between the wavevector com- 
ponents k x , k y and the frequency u>. Below the plasma 
frequency, the permittivities of metal and dielectric have 
different signs, and this ensures a key feature of metal- 
dielectric nanostructures (MDNs): the principal elements 
of their permittivity tensor can have nearly arbitrary val- 
ues. For example, if £\ — —e 2 d\jd 2 then the structure 
has very large permittivity (ej_ — > 00), and it repre- 
sents realization of an "cpsilon-vcry-large" (EVL) ma- 
terial [13]. If £1 = —e 2 d 2 jd\ then an MDN has near 
zero permittivity (£|| = 0) and is a good candidate for re- 
alization of the so-called "cpsilon-near-zero" (ENZ) ma- 
terial [14]. Basically, the EVL and ENZ materials are 
optical conductors and insulators which form a basis of 
metactronics [10]. If £y and e± have different signs then 
the MDN is a typical realization of the so-called indefinite 
medium [15]. 

In this work we considered three MDNs formed by lay- 
ers of metal and dielectric with various thickness ratios, 
but fixed total period. Configurations of the structures 
are illustrated schematically in the insets of Figs. 2(a-c). 



The dispersion diagrams D/X(k y ) for the three MDNs 
under consideration shown in Figs. 2(a-c) were computed 
using two approaches: the effective medium model (ap- 
proximate approach) and the well-known classical disper- 
sion relation for one-dimensional photonic crystals (exact 
description). 

The properties of such MDNs are investigated at wave- 
lengths which are approximately 10 times larger than the 
structure period. The permittivity of dielectric is as- 
sumed to be constant and equal to e\ = 4.6. Permittivity 
of metal is given by the Drude model: e 2 = 1 — A 2 /Ap, 
where X p = 4D. The losses in metal are neglected in 
order to simplify the consideration. The frequency de- 
pendencies of £|| and e± for three types of MDNs with 
different ratios of layer thicknesses are shown in Fig. 1. 
The singularities of £j_ and nulls of £|| corresponding to 
EVL and ENZ behavior, respectively, are clearly visible. 

The information about dispersion properties of MDN 
provided by effective medium model is approximate, and 
it does not take into account actual periodicity of the 
structure. The general expressions for s« and e±_ in 
Eq.(l) are deduced under quasi-static approximation by 
using an assumption that the electric field does not vary 
inside of the layers. If permittivities of the constituent 
layers are positive then the requirement that the period 
of the structure (and thus the thicknesses of the con- 
stituent layers) is much smaller than the wavelength of 
operation ensures negligibly small variation of fields since 
spatial harmonics propagating inside of the layers have 
small wave numbers. However, such a reasoning is not 
valid if materials with negative permittivity (e.g. met- 
als) are involved because the interfaces between materi- 
als with positive and negative permittivity support sur- 
face plasmon polaritons (SPPs) which form an alterna- 
tive channel for wave propagation and may cause sig- 
nificant field variation in the layers. The SPPs at the 
individual interfaces of an MDN exist at the frequencies 
below the SPP resonance which corresponds to condition 
£1 = —£2. It is logical to expect violation of classical ho- 
mogenization requirement at these frequencies, but only 
for oblique propagation of waves through the MDN since 
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the waves propagating across the layers cannot couple to 
SPPs. Note, that the SPP resonance condition, which 
corresponds to the strongest confinement of SPPs and 
thus leads to strongest variation of fields across the lay- 
ers, connects only permittivities of the materials which 
form the interfaces where SPPs propagate, and it does 
not depend on the thicknesses of the layers. 



The dispersion characteristics of actual periodic MDN 
can be evaluated from Maxwell's equations with imposed 
periodic boundary conditions. A consideration of the 
MDN as one-dimensional photonic crystal leads to the 
following exact dispersion relation for eigenmodes with 
transverse magnetic (TM) polarization: 
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cos(k x D) = cas&M) cos(fc( 2 )d 2 ) - - ™ + sinCfcW^) sin(fc( 2 >d 2 ), (3) 
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where k£ = [zi,2{u / c) 2 — k 2 ,] 1 / 2 are x-components of wavevectors in the corresponding layers. 



The MDNs with 3 different ratios of layers thicknesses 
(3/2, 1 and 2/3, respectively) were chosen in order to 
demonstrate different behaviors of dispersion curves. In 
the first case, the SPP resonance (in our case correspond- 
ing to D/X = 0.1055) lays over the frequency where 
EVL behavior is expected {D/X = 0.089). In the sec- 
ond case, the frequencies are chosen to be equal. In the 
third case, the SPP resonance appears well below EVL 
frequency {D/X = 0.124). In all cases the dispersion 
curves consists of two branches and have SPP resonance 
as asymptote if k y — > oo. However, in the 1st case the 
upper branch corresponds to a backward wave while the 
second branch doesn't. In the 2nd case both branches 
correspond to forward waves and co-exist at the same 
frequency range in contrary to the 1st case. In the 3rd 
case the branches cross each other right at the frequency 
corresponding to ENZ behavior of MDN {D/X = 0.089) 
and one of the branches has a maximum leading to ex- 
istence of backward and forward waves simultaneously 
at the same range of frequencies. The effective medium 
model in all cases predicts only one propagating wave at 
all frequencies. The presence of two propagating waves 
is a consequence of nonlocality and strong spatial disper- 
sion which are caused by SPPs at the interfaces of the 
layers. 

We have verified that the behavior observed for MDNs 
with the period 10 times smaller than wavelength re- 
mains qualitatively the same for MDNs with 2, 4 and 
8 times smaller periods (i.e. 20, 40 and 80 times smaller 
than wavelength) . The convergence to an asymptotic line 
corresponding to SPP resonance becomes weaker (the 
branches reaches this line at larger k y ) , but the topology 
of the branches remains the same. On the other hand, the 
decrease of period pushes the dispersion curves closer to 
those predicted by the effective-medium approximations 
for small k y . It is interesting that this convergence also 
depends on the ratio of the layers thicknesses. In the case 
of equal layers the convergence is faster than for different 
layers. However, even for very small period the corre- 




FIG. 3: (Color online) Diagrams of beam refraction at the air- 
MDN interface plotted both for local effective medium model 
and actual MDN with help of isofrequency contours at three 
different frequencies: a) D/X = 0.085, (bj D/X = 0.094, and 
(c) D/X = 0.106. 

spondence between the results of the effective model and 
actual dispersion properties of MDNs is not satisfactory, 
especially for the waves with relatively large k y which are 
mainly used in the problem of subwavclcngth imaging. 

The isofrequency contours can be used for describing 
the beam refraction at the interfaces separating free space 
and the layered structure. They allow to understand the 
character of refracted waves and find directions of their 
phase and group velocities. For the three typical cases of 
isofrequency contours, the refraction diagrams are plot- 
ted in Fig. 3. In the first case [sec Fig. 3(a)], there exists 
only one negatively refracted beam, both for the effec- 
tive medium model and transfer-matrix method. Here 
we have a situation where the effective medium model 
can be applicable. The same situation is demonstrated in 
Fig. 3(c), but one positively refracted beam is observed. 
The simultaneous appearance of ellipse and hyperbolic- 
like contour in Fig. 3(b) leads to birefringence phenom- 
ena in the metal-dielectric structure. In this case both 
negatively and positively refracted TM-polarizcd beams 





should appear simultaneously. 

In order to verify this prediction we have performed 
numerical simulations with a commercial full-wave elec- 
tromagnetic solvers package [16] for the wave refraction 
at an interface at different frequencies [sec Figs. 4(a-c)]. 
It should be noted that such simulations are quite hard 
to accomplish since the beam width is equal to about 
10 wavelengths of radiation A, with the smallest mesh 
element reaching A/500. At D/X = 0.085 only one neg- 
atively refracted wave is observed [Fig. 4(a)] and the ef- 
fective medium model can be applicable. Appearance of 
an ellipse in the isofrequency contours at D/X = 0.094 
shown in Fig. 3(b) leads to birefringence [Fig. 4(b)]. Neg- 
atively refracted beam still appears but in addition to it 
a positively refracted beam is observed. In this case, the 
effective medium model is not applicable, since it does 
not predict the presence of the two waves describing only 
one of them. At the higher frequency (D/X — 0.106) only 
one positively refracted ray is observed, and this fact is 
well described by the effective medium model. 

Our results call for deeper studies of the applicabil- 
ity limits of the effective medium models, and they may 
be useful to explain experimental results in more complex 
settings. In particular, a few years ago Kozyrev et al. [17] 
observed experimentally the waveform splitting into two 
or multiple beams after the propagation of electromag- 
netic waves through a slab of a magnetic metamaterial. 
Such effects can not be described by an effective-medium 
model being characterized by some effective permittiv- 
ity and permeability. Instead, the observed transmission 
properties of metamaterials are affected significantly by 
the internal structure, such that the multiple beam for- 
mation observed in the experiments can be attributed to 
the excitation of magnetoinductive waves. 

In conclusion, we have studied the dispersion prop- 
erties of periodic metal-dielectric structures by employ- 
ing the exact transfer-matrix approach and an effective 
medium model. We have revealed a substantial difference 
between the results provided by these two approaches, 



even in the limit when the structure spacing is much 
smaller than the radiation wavelength. These result ques- 
tion the applicability of the averaged parameters and ef- 
fective media in the cases when the structure support sur- 
face waves. In particular, we have found two dispersion 
branches of extraordinary waves instead of one branch 
predicted by the effective medium model. Importantly, 
the effective optical nonlocality has been shown to de- 
pend dramatically on a ratio between the thicknesses of 
metal and dielectric layers, and it can be engineered by 
changing this ratio, with the weakest nonlocality in the 
case of equal thicknesses. Strong nonlocal effects pre- 
dicted here can be observed as the splitting of the TM- 
polarized wave at an interface between air and metal- 
dielectric nanostructured metamaterial into two waves 
with positive and negative refraction. Our results sug- 
gest the presence of strong spatial dispersion in many 
types of nanostructured materials and, in particular, they 
provide a proof that metal-dielectric nanostructures are 
inherently nonlocal metamaterials. 
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